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I. INTRODUCTION 



We report on the results of a calculation of the next-to-leading (NLO) order QCD cor- 
rections to the QCD processes 

pp -> W ± +2 jets -> £ ± u +2 jets , 

pp -> Z/ 7 * + 2 jets -> + 2 jets . (1) 

(single flavor, I = e, /x or r) at the CERN Large Hadron Collider (LHC), i.e. pp collisions 
at y/s = 14 TeV. We ignore electroweak (EW) contributions arising from e.g. t-channel 
W/Z/j exchange between two scattered quarks, which are typically 0(1%) for total rates. 
Calculations for the Fermilab Tevatron, pp collisions at y/s = 2 TeV, were considered in a 
previous publication [1], which also introduced the calculational technique. The results are 
made with the Monte Carlo program MCFM, which makes full predictions for any infra-red 
safe variable, including fully differential distributions, for any set of experimental cuts, and 
any decay modes of the Z/j* intermediate states. However, here we consider only leptonic 
vector boson decays. All leptons are treated as massless, so the results treat W, Z — > e, //, r 
flavors as equivalent. 

MCFM uses previously published matrix elements for the crossed reactions e + e~ — > 
4 partons [2] and e + e~ — > 5 partons [3], in the four dimensional helicity scheme. The real 
corrections to the basic Born processes, i.e. 

parton + parton — > W/Z/^f* + 3 partons (2) 

were published in Refs. [3-5]. MCFM incorporates these matrix elements using the subtrac- 
tion method [6, 7]. These details and those of the numerical checks performed may be found 
in Ref. [1]. 

As a useful comparison to more inclusive, well-studied weak boson production, we simul- 
taneously present MCFM results for the processes W, Z + 0,1 jets. QCD corrections for these 
processes have been known for some time [8], and have already proved invaluable for Teva- 
tron studies. The inclusion of the W/Z + 1 jet processes in MCFM aided in understanding 
the issues to be faced in implementing the more complicated W/Z + 2 jet processes. The 
next-to-next-to-leading order (NNLO) corrections to total inclusive W, Z production are 
known [9], but are not implemented in MCFM. 
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Parameter 


Default value 


Parameter 


Default value 


M z 


91.187 GeV 


a(M z ) 


1/128.89 


T z 


2.49 GeV 


G F 


1.16639xl0~ 5 


M w 


80.41 GeV 


al 


0.42662 (calculated) 


T w 


2.06 GeV 


sin 2 9 W 


0.23012 (calculated) 



TABLE I: Default parameters in the program MCFM. 

We also present results where the additional jets are a heavy flavor pair, specifically b 
quarks, i.e. W ± bb ) Zbb, and compare to the flavor-inclusive expectations. This was previ- 
ously studied for the Tevatron [10, 11], but our results are the first NLO predictions for 
these rates at the LHC. These are of particular importance as backgrounds to a number of 
new physics searches, especially for the Higgs boson. Experimental studies have so far had 
to rely on leading order (LO) predictions. For the W ± bb processes we find a LO calculation 
underestimates the rates at LHC energies by about a factor 2.4. 

II. GENERAL RESULTS 

Here we present results for the total observable cross sections for the production of a 
vector boson and up to 2 jets. Results specific to Higgs boson searches can be found in 
Section IV. 

A. Input parameters 

For all our results we use the default set of EW parameters in MCFM, which are given 
in Table I. As noted in the table, some parameters are calculated using the effective field 
theory approach [12], 

G e 2 
e 2 = A7ra(M z ), gl = 8M^-^, sm 2 6 w = —. (3) 

For simplicity we have taken the CKM matrix to be diagonal in the W + 2 jets process. As 
a consequence there are, for example, no annihilating us initial states for this case. This 
approximation is not expected to influence any anticipated analyses. For the other processes 
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we retain only the Cabibbo sector of the CKM matrix, 



( 0.975 0.222 0^ 



V, 



CKM 



0.222 0.975 
1 



(4) 



The value of as(Mz) is not adjustable; it is determined by the chosen parton distribution. 
A collection of modern parton distribution functions (PDFs) is included with MCFM, but 
here we concentrate on the CTEQ6L1 (LO) and CTEQ6M (NLO) sets [13] which specify 
a s (M z ) = 0.130 and a s (M z ) = 0.118 respectively. Unless otherwise stated, we take the 
factorization and renormalization scales to be ji = /if = /i r = My. 



B. Basic kinematic cuts and jet selection 

For the general results presented here, we consider single leptonic decays, 

W + -> ve + , W~ -> ue~ , Z/ 7 * -> e - e + . (5) 

W + and W~ must be considered separately for pp collisions, as their rates are not equal, 
given the unequal quark/anti-quark parton distributions. 

Since these results are also intended to help calibrate experimentalists' Monte Carlo 
results, the results must be obtained with kinematic cuts on the final state particles for 
which they will be observable at the LHC experiments ATLAS and CMS. All leptons must 
satisfy 

p T (£) > 15 GeV, \ m \ < 2.4, (6) 

We also require that the charged dilepton invariant mass be greater than 15 GeV. This 
prevents the production of soft e~e + pairs which would otherwise be copiously produced 
by the virtual photon in the Z/Y process. We use the Run II kr clustering algorithm [14] 
to find jets, with a pseudo-cone of size R = 0.4. Jets are also subject to kinematic cuts 
consistent with detector requirements for observability: 

PtO') > 20 GeV, 1^1 < 4.5. (7) 

Finally, we require a minimum separation to isolate the leptons: 

AR Ej > 0.4 , AR U > 0.2 . (8) 
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For the specific case of heavy-flavor jets, which are experimentally observed via vertex tag- 
ging, the pseudo-rapidity restriction is more severe: 

\Vb\ < 2.5 . (9) 

We do not impose a cut on missing transverse momentum, f Tl for the W ± + jets cases. Not 
all analyses planned for the LHC use such a cut, or use different values for it, depending on 
what else is required in the event and machine luminosity. Instead, we show the differential 
distribution in this variable. 

Except where noted, we consider only inclusive jets production: W, Z + n jets refers to 
the cross section for production of W, Z and n or more jets. 

C. Scale dependence 

The principal motivation for performing a NLO calculation is to reduce the uncertain- 
ties in LO predictions. In particular, any perturbative prediction contains an unphysical 
dependence on renormalization and factorization scales fJ> r ,fJ>f. The scale /if is introduced 
during the factorization of the calculation into a perturbative hard scattering part and 
non-perturbative PDFs. The latter are taken as input from data, with additional perturba- 
tive evolution. The factorization and renormalization scales are often chosen to be equal, 
/i = /i r = /if, and we do this here. Strictly speaking, the two scales are unrelated and a 
calculation should be checked for independence of each scale individually. However, MCFM 
does not currently allow for this possibility for most processes. 

The magnitude of cross sections and the shape of differential distributions can vary greatly 
between two different choices of scale. This is often interpreted as an inherent "theoretical 
uncertainty". Another strategy is to argue for a particular choice of scale, based on the 
physics of the process under consideration. Sometimes this gives results at LO that are 
close to the NLO results, but is not guaranteed. 

A NLO calculation is an invaluable tool for investigating the issue of scale dependence. 
In a calculation performed through to order a™, the residual scale dependence enters only at 
order a™ +1 . As a result, one expects that NLO predictions are more stable under variations 
of the scale, \i. In addition, the NLO result may provide further evidence to support a 
particular scale choice that may have been deemed appropriate at leading order. In this 
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FIG. 1: Scale dependence of the eu e + jets cross sections, fi = \i T = fif, using the cuts described in 
Sec. II B. In general, ^-dependence is significantly reduced at NLO compared to LO. 

case, the NLO result does not imply that the scale choice is physically meaningful, but it 
can provide an easy method of normalizing a LO calculation. (A LO calculation typically 
proceeds much more rapidly, and may need to be repeated multiple times for different 
kinematic configurations.) 

To show the improvement achieved in the present calculation, we plot the scale depen- 
dence of the total inclusive W, Z + 0, 1, 2 jets cross sections, including leptonic decay, with 
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FIG. 2: Scale dependence of the e + e~j and e + e~jj cross sections, /x = ix r = fif, using the cuts 
described in Sec. II B. 

the kinematic cuts of Eqs. (6-8), in Figs. 1,2. The NLO predictions for W, Z + 0, 1 jets 
have been known for some time [8, 15, 16], but we recalculate them here with MCFM. The 
W + 2 jets behavior, Fig. 1, is practically independent of //, as in the W, Z + 0,1 jets cases, 
demonstrating a small theoretical uncertainty. Table II summarizes the NLO results for 
V+ jets processes, along with the residual theoretical uncertainty due to scale dependence. 
All Monte Carlo statistical uncertainties are less than 1%. While the QCD corrections for 
V + 0,1 jets are positive, for V + 2 jets they are slightly negative. 



D. Kinematic distributions 



Once again, we repeat some W, Z + 1 jet results, in order to highlight both the similarities 
and the differences with the corresponding 2-jet distributions. We first show the missing 
transverse momentum distributions for VF 1 *^ jets in Fig. 3, and then the leading (hardest) 
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"z/ e + X 
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3970 


4830+S 

— yu 


^ 


- e ~ +X 


803 


915 ±31 




'v e j + X 


1660 


18801™ 


e" 


~v e j + X 


1220 


1420 ± 40 


e H 


-e-j + X 


248 


288+^ 


e H 


~v e jj + X 


773 


669l? 8 


e" 


~&ejj + X 


558 


4911? 


e H 


-e~ jj + X 


116 


105±S 



TABLE II: Summary of LO and NLO cross sections [pb] for W/Z + 0, 1,2 jets, including single 
leptonic decay of the weak boson. The central value at NLO is for fi = My, and the uncertainty 
at NLO is for scale variation from \i = My/2 to \i = 2My. No uncertainties are shown for LO 
results, as they are all considerably larger than at NLO. All Monte Carlo statistical uncertainties 
are less than 1%. 

jet p T distribution for the W ± j and W ± jj cases in Fig. 4, and for Zj and Zjj in Fig. 5. 
The leading order curves for W ± + 1 jet in Fig. 3 show a kinematic structure not present in 
the other curves. At leading order, the cut on the minimum px of a jet, Eq. (7) also imposes a 
minimum transverse momentum for the W boson, since at this order the transverse momenta 
are equal and opposite. Most events lie right above the cut, Pt{W) ~ 20 GeV. The vector 
addition of this peaked distribution with the transverse momentum generated in the decay 
leads to the structure shown in Fig. 3. At NLO (or for higher jet multiplicity), the jet 
transverse momentum cut does not impose the same constraint on the momentum of the W 
boson. This explains the absence of the structure in the other plots. The existence of this 
structure in the LO plots reinforces the importance of NLO calculations. 

In contrast to the W, Z + 1 jet cases, the px distributions turn over at small values for 
V + 2 jets. If the highest pt jet has a transverse momentum, p^,\ close to 20 GeV, there is 
little phase space for the emission of a second, softer jet with p^ > pjp > p™ m = 20 GeV. 

We show the boson pt distributions in Figs. 6,7, and the dijet invariant mass distributions 
in Figs. 8,9. The former show a pronounced NLO deficit relative to LO at large values of 
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Pt{V), whereas for the latter there is almost no change in shape from LO to NLO. The same 
holds true for the pseudo-rapidity separation of the jet pair, as shown in Figs. 10, 11. 



III. HEAVY FLAVOR CONTENT OF JETS 

We calculate the fraction of W, Z + 2 jet events that contain two heavy quark jets, limiting 
ourselves to b quarks, because they can be tagged with high efficiency in experiment. These 
processes are of particular interest as backgrounds to new physics searches, in particular 
Higgs bosons. The calculations are identical to those presented in Ref. [1], but for pp 
collisions at the LHC energy of ^/s = 14 TeV. We work in the approximation in which the 
b quarks are taken to be massless, as the massive results are not yet known at NLO, and 
we ignore contributions from processes in which there are b quarks already present in the 
initial state. 
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FIG. 4: Hardest jet pr distribution of e ± u e j and e ± u e jj events at LO (dashed) and NLO (solid). 
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FIG. 5: Hardest jet pt distribution in e + e~ j and e + e~ jj events at LO (dashed) and NLO (solid). 

We begin as in the non-heavy flavor case, checking the theoretical uncertainty of the 
calculation by plotting the scale dependence of the cross sections, shown in Figs. 12,13. 
The W + bb cross section shows much greater \x dependence at NLO than at LO, contrary to 
naive expectations. For the central scale choice \x = M W) the NLO result is a factor 2.4 larger 
than at LO, and requires some explanation. At LO, charge conservation allows only qq initial 
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FIG. 12: Scale dependence of the e^v e bb cross sections, /x = n T = [if, using the cuts described in 
Sec. II B. Note the behavior of Wbb at NLO: the NLO rate with initial gluon density set to zero 
is well-behaved, but with real initial gluons it behaves much more like a LO calculation. See text 
for additional comments. 



Oh 



pp^e e bb+X 



NLO 



LO 



50 



100 

fi [GeV] 



150 



200 



FIG. 13: Scale dependence of the e + e~ bb cross section, jjl = fj, r = fj,f, using the cuts described in 
Sec. II B. The Zbb result is better behaved than in the case of Wbb, but still retains an uncertainty 
of order ±15%. 

states to contribute to Wbb production. At NLO, the same is true for the virtual corrections, 
but not for the real emission corrections, which contain subprocesses qg — > q'Wbb. Some 
representative Feynman diagrams are shown in Fig. 14. This allows the very large gluon 
luminosity at the LHC to contribute. However, if we ignore the additional bb pair, then NLO 
corrections are of Drell-Yan type, Fig. 14(a), and might be expected to receive an additional 
contribution of only 30% at LHC energies (see Table II). The remaining contribution is 
due to diagrams of the type shown in Figs. 14(b,c), which contain no singularities for a 



13 



Diagrams by MadGraph 






(a) (b) (c) 

FIG. 14: Some representative diagrams for the process ug — > dW + bb. See text for discussion. 
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component. 

finite b mass; W emission requires the massless quark propagator to be off mass shell [17]. 
(We calculate with m& = 0, but avoid the massless singularity by imposing a finite cut 
on both b quarks.) The contributions from diagrams of the type in Fig. 14(b) open up 
additional phase space, since they are not suppressed by an s-channel gluon propagator 
which occurs in diagrams of the type shown in Fig. 14(a). Furthermore, contributions from 
diagrams of the type in Fig. 14(c) contain a large Casimir factor from the three-gluon vertex. 
Diagrams (b) and (c) constitute new, LO contributions that cannot appear in Drell-Yan. 
Note, however, that one cannot consider individual diagrams separately, as they must all be 
included together at the amplitude level to maintain gauge invariance. We demonstrate that 
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the additional diagrams are the cause of the hugely enhanced cross section by calculating 
the NLO W ± bb cross section with initial gluons turned off, i.e. g(x) = 0. This result has 
practically no \i dependence, as shown by the dotted curves in Fig. 12. This is further 
supported by examining the rapidity distribution of the emitted W boson, shown in Fig. 15. 
The significant NLO enhancement at large rapidity, resulting in a noticeable plateau in the 
distribution, demonstrates the large size of the gg-initiated component: quarks on average 
occur at larger values of Feynman x than gluons, resulting in an overall boost of the W 
boson in the rest frame of the detector. Thus, the dominant component of inclusive W^bb 
production at LHC energies, the summed gluon-initiated real emission subprocesses, still 
contains large LO theoretical uncertainty. This can be resolved only by performing a NLO 
calculation of W ± bbj. 

This situation does not exist for Zbb production, because lack of a flavor- changing con- 
straint on any quark line means that initial state gluons contribute already at LO. This is 
readily apparent from Fig. 13, which shows NLO corrections much closer to the size expected 
for Drell-Yan processes, and a residual uncertainty due to scale dependence of about 20%. 
Note the significant feature of the heavy flavor case, that a(W + bb) > a(Zbb) > a(W~bb), 
as opposed to <i(W + jj) > o , (W~jj) ^> a(Zjj) for the flavor-inclusive cross sections. 

As for the W ± jj case in Sec. II D, we histogram the differential W^bb cross sections as a 
function of the p T of the leading jet, shown in Fig. 16, along with the desired ratios of W ± jj 
to W ± bb cross sections, to illustrate the relative behavior of heavy flavor production. In 
principle, W ± cc production would behave similarly to W ± bb. We also show the distributions 
as a function of dijet invariant mass in Fig. 17. Including the b quark mass at LO is only 
a 4% effect on the total cross section, but is contained almost completely in the invariant 
mass range 2rrib < rribb < 20 GeV. Consequently, the two curves would exactly overlap 
each other in Fig. 17, except in the first bin. 

Similar histograms for the Zbb case are shown in Figs. 18,19. Table III summarizes the 
W ± , Z + bb cross section results and scale uncertainties. 

IV. APPLICATIONS 

One of the principle goals of the LHC program is to determine the nature of electroweak 
symmetry breaking (EWSB), which in the Standard Model (SM) is due to a complex scalar 
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TABLE III: Summary of LO and NLO cross sections [pb], including the theoretical (scale) uncer- 
tainty at NLO, for W/Z + bb, including leptonic decay of the weak boson. The central value at 
NLO is for fi = My. All Monte Carlo statistical uncertainties are less than 1%. 

577(2) doublet field. Upon acquiring a vacuum expectation value, it gives rise to both 
weak boson and SM fermion masses and creates a single physical scalar, the Higgs boson. 
While there are other variations on this idea, or other models which generate EWSB, nearly 
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FIG. 17: The dijet invariant mass distribution of e^v e bb events, and of the ratio of W^jj to 
bb, at LO (dashed) and NLO (solid). For the latter, the jets from W ± jj are also restricted to 
the region \yj\ < 2.5, to match the acceptance cuts that will be used for bb events. 
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FIG. 18: The leading jet pr distribution of e + e~ bb events, and of the ratio of Zjj to Zbb, at 
LO (dashed) and NLO (solid). For the latter, the jets from Zjj are also restricted to the region 
\yj\ < 2.5, to match the acceptance cuts that will be used for Zbb events. 
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FIG. 19: The dijet invariant mass distribution of e + e~ bb events, and of the ratio of Zjj to Zbb, 
at LO (dashed) and NLO (solid). For the latter, the jets from Zjj are also restricted to the region 
\yj\ < 2.5, to match the acceptance cuts that will be used for Zbb events. 

all of them have a physical scalar field which couples to weak bosons and SM fermions 
similar to the SM Higgs boson. Many searches for it at hadron colliders look for associated 
Higgs boson production, pp — > W ± H,ZH [19], or a Higgs boson produced between two 
far forward/backward high-p^ jets. The latter is weak boson fusion (WBF) Higgs boson 
production, which arises from a pair of initial-state quarks both emitting either a W + W~ 
or ZZ pair, which fuse to create the Higgs boson [20]. 

Precise theoretical knowledge of the W, Z+ jets cross sections is of particular interest 
at the LHC, since these processes constitute significant backgrounds for the above Higgs 
boson searches, as well as other SM processes of interest, such as top quark production, and 
many other new physics searches. We choose to examine only two particular cases: W^bb 
as a background to W ± H; H — > bb, and Zjj as a background to WBF Higgs production, 
pp -> Hjj. 



A. W±bb and W ± H: H —> bb 

The most desirable Higgs boson search channels at the LHC are inclusive production 
(gg — > H via a top quark loop) or WBF Higgs production, due to their much larger cross 
section and much cleaner background environment, respectively. Despite the smaller rate, 
the W ± H; H — > bb channel has generated great interest because of its potential to measure 
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FIG. 20: The 66 invariant mass distribution of e ± f e 66 events at LO (dotted), NLO (dashed), and 
NLO with a jet veto on the non-6 real emission (solid), with kinematic cuts specific to a potential 
W ± H; H — > bb search as described in the text. 

the bottom quark Yukawa coupling. Several LHC experimental groups [21, 22] have studied 
the channel in detail. These studies claim that the search will be quite difficult for Higgs 
boson masses above the current LEP limits, but possible. While the signal is known at 
NLO, our calculation of the W^bb background at NLO is the first to rectify this omission. 

As shown in Sec. Ill, the NLO cross section is more than two times the LO result. This 
does not bode well for this particular Higgs boson search channel. However, we must apply 
the same kinematic cuts as in the current analyses and examine the NLO/LO ratio in a dijet 
invariant mass bin relevant for the search. We apply the following cuts as in Ref. [21]: 

p T (£) > 20 GeV , \y e \ < 2.4 , (10) 
p T (j) > 30 GeV , |%|<2.5, (11) 
AR j£ > 0.4 , AR U > 0.2 , (12) 

where the jet may or may not be a b jet. The resulting bb invariant mass distribution is 
shown for the summed W + and W~ components in Fig. 20. Assuming a Higgs boson mass 
of Mjj = 120 GeV, and a mass bin of ±20 GeV, i.e. 100 < < 140 GeV, the NLO cross 
section is a factor 1.9 larger than LO. 

Since we know that the component causing the large enhancement typically has an ad- 
ditional jet emission at large pr, we impose a non-6 jet veto at NLO in the real emission 
calculation. This turns the cross section into an exclusive result, where the heavy flavor is 
assumed to have been tagged. This NLO exclusive result is about 10% smaller than the LO 
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result, making the background much more manageable. However, we emphasize that use 
of this result requires normalizing a parton shower calculation to this cross section after a 
light-jet veto is imposed in the parton shower Monte Carlo! The residual uncertainty from 
scale variation is ~tyo %• 

B. Zjj and Hjj 

A characteristic feature of WBF Higgs boson production is Hjj events with the Higgs 
boson appearing at central rapidities, but the two scattered quarks appearing at large rapidi- 
ties with significant pr- Because the process is color singlet exchange, the events typically 
have little central jet emission. In contrast, QCD Zjj production prefers the jets to be much 
more central, while the Z comes from Bremsstrahlung and is at larger rapidites. Also, the 
color triplet or octet exchange gives rise to significant additional central jet emission, known 
as "minijet" activity [23, 24]. We use MCFM here to provide a NLO normalization of the 
QCD Zjj rate in a region of phase space relevant for WBF Higgs boson searches where the 
Higgs boson decays to tau pairs. 

Standard kinematic cuts for these searches will require the jets, frequently called "tagging 
jets" , to have a minimum separation of about 4 units of pseudo-rapidity, to appear in opposite 
hemispheres, and for the Higgs decay products (e.g. t + t~) to appear between the jets [25]. 
In addition to the basic cuts of Eqs. (6-8), we thus impose the following cuts: 

\r,f ] - r,f\ > 4.0 , r,f ■ r,f ] < , rjf n + 0.4 < m , M < V ™ x - 0.4 , (13) 

where r^ ax (r]j im ) is the larger (smaller) of {v^ivf^}- This reduces the QCD Zjj cross 
section from 94.9 pb to 2.9 pb at LO, and 104 pb to 3.4 pb at NLO. Comparison to previous 
experimental studies are difficult, as all of these used dynamical scales, and were all different 
from each other [26]. 

We show the scale dependence for Z + 2 jets in this phase space configuration in Fig. 21. 
The NLO is quite stable, as in the total cross section result, except for very low values of \x. 
This is not uncommon in NLO calculations and typically indicates the presence of logarithms 
that can become large with small choices of fi, but are not necessarily a problem at that 
order in the perturbative expansion. We also plot several differential distributions in Fig. 22 
to demonstrate that there do not appear to be any large NLO corrections to the shapes in 
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FIG. 21: Scale dependence of the e + e~jj (t + t~jj) cross section, /x = ^ r = /xy, using the cuts 
described in Sec. IV B to examine the region of far forward tagging jets. 



this small region of phase space; the NLO distributions in general appear to be very similar 
to the LO distributions. One notable exception is the azimuthal angle between the tagging 
jets, which is flatter at NLO, compared to the preferential back-to-back configuration of the 
jets at LO. This is likely due to emission of final state radiation spreading out the preference 
for the QCD jet pair to be back-to-back, but requires further investigation. This distribution 
is of special relevance for studies of Higgs CP transformation properties [27]. 



V. CONCLUSIONS 

We have presented the first results for the implementation of W/Z + 2 jet production 
at NLO at the LHC, including the specific cases of heavy flavor jets. An analysis based on 
inclusive jet production for the LHC shows that in most cases the usual benefits of NLO are 
realized, among them a reduced scale dependence and hence an improved normalization for 
kinematic distributions. The outstanding exception is the case of W ± bb production, which 
at one higher order in the QCD coupling experiences additional LO contributions that 
dramatically enhance the total rate. We showed that this additional contribution typically 
leads to an additional non-6 jet at large pt and pseudo-rapidity, which can be identified in 
the detectors. 

These processes are of particular interest at the LHC as backgrounds to other SM pro- 
cesses and searches for new physics. We analyzed two cases in particular. The first analysis, 
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FIG. 22: Various kinematic distributions for QCD e + e jj (r + r jj) production as a background 
to Hjj production, as described in the text. 
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W ± bb as a background to WH; H — > 66 for a light Higgs boson, suffers considerably from the 
unexpected large NLO corrections, but may be controlled somewhat by vetoing non-6 jets. 
This will require additional analysis with parton shower Monte Carlo normalized to our NLO 
jet-vetoed rates, after the veto of additional jets radiated in the parton shower. The second 
analysis, Zjj as a background to weak boson fusion Higgs production, Hjj, shows that the 
overall corrections are not very large for a central scale choice of \i = Mz, only about +14%. 
However, most Hjj analyses have used dynamical scale choices of pr(j, min), which average 
only a few tens of GeV, so these analyses will generally benefit from a reduction in the QCD 
Zjj background normalization. In addition, we have shown that no large NLO corrections 
appear in this tiny region of phase space with jets widely separated in pseudo-rapidity, and 
various kinematic differential distributions at NLO appear very much like their LO counter- 
parts. This allows simple overall normalization of parton shower Monte Carlo based on the 
total rates, but does not address the issue of vetoing events with additional minijet activity. 
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